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Abstract. We introduce some analytic relations on the set of partial differ- 
ential equations of two variables. It relies on a new comparison method to 
give rough asymptotic estimates for solutions which obey different partial dif- 
ferential equations. It uses a kind of scale transform called tropical geometry 
which connects automata with real rational dynamics. Two different solutions 
can be considered when their defining equations are transformed to the same 
automata at infinity. We have a systematic way to construct related pairs of 
different partial differential equations, and also construct some unrelated pairs 
concretely. These verify that the new relations are non trivial. We also make 
numerical calculations and compare the results for both related and unrelated 
pairs of PDEs. 



1. Introduction 

Scaling limits play one of the central roles in discrete dynamical systems, which 
create another dynamics with different analytic aspects. Tropical geometry and 
ultradisctete dynamical systems arose from very different contexts each other, and 
have been developed quite separately [V], [TTMS]. However from the view point 
of scaling limits, their method surprisingly coincides, which connect real rational 
dynamics with piece-wise linear systems. 

Such scaling limit provides with an important prototype in mathematical physics. 
The Korteweg-de Vries (KdV) equation Ug + uux + u^x = is a PDE which has 
been derived from the discrete Lotka-Volterra equation z*'^^^ = + -|- 

■^Af+2)~^ [Hi IHT] . The above scaling limit changes the discrete dynamics into the 
Lotka-Volterra cell automaton t/^^ = ^Ar+i + inax(0, U]^^) — max(0, In 
the light of integrable systems, these three dynamics share common features, like 
existence of soliton solutions, many conserved quantities and so on [TTMS]. 

These scaling limits allow us to study analytic aspects at the same time for 
three categories of dynamical systems which sit in different hierarchies mutually. 
The classes of dynamics which can be analyzed by these scaling limits are rather 
broad beyond integrable systems. In fact it is easy to find discrete dynamical 
systems which are far from integrable, but which are transformed into integrable 
cell automata, (see section [SIC). If we consider a situation when two discrete 
dynamical systems are transformed into the same integrable cell automaton, where 
one is integrable and the other is not, then it will be quite natural to try to produce 
some systematic method to analyze various PDEs which include wide classes of 
dynamical systems as above. 
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Arising from a new comparison method to study rough asymptotic growth among 
solutions to different partial differential equations, in this paper we introduce ana- 
lytic relations on the set of partial differential equations of two variables. As a first 
step in this paper we show non triviality of such relations. 

Let -Pi(e, u,Ux,Us, ■ ■ ■) and P2{e,v,Vx,Vs, ■ ■ ■) be two different partial differential 
equations for dependent variables u{x, s) and v{x, s) respectively. Our basic method 
is to compare solutions between Pi and P2 in terms of the initial data and the 
distorsions of the higher derivatives of their solutions. 

For small parametrization e > 0, let: 

PDE2 = {P(e, u, Ux, Us,Uxs, • ■ • , Uas) = : 

P are polynomials with real coefficients} 

be all the set of families of polynomial type partial differential equations with 2 
variables parametrized by e > 0. 

In this paper we introduce some analytic relations ^ between elements in PDE2. 
They are given by uniform asymptotic estimates for all positive solutions of 
different PDEs. We will focus on two aspects, where one is growth rates of solutions 
and the other is the domains of solutions on the space variable. For the first we will 
treat two growth cases, exponential and double-exponential . For the second 
we also treat two cases, where one is which include domains of bounded space 
variables, and the other is '^oo which restricts only on infinite domains in space. 
Thus in total there are four classes of the relations: 

fin — fin 

|A |A 



where ^+ < -^^z means that u ^+ v implies u v for any two solutions. 

In this paper we introduce two subsets: 

PDE^ " C PDE^ C PDE2. 

There are also some stratifications (M, c, D, L; a) and (M, c, L, k, D) over both ~+ 
and PDE2 with respect to some constants respectively (2.B,3.B). 

Our main aim here is to verify non triviality of the relations over these spaces. 
In [K2], we have introduced some constructive way to obtain partial differential 
equations, which produces many related pairs with respect to the classes: 

(pde/-,~},j, (pde-,-:,) 

where * = e or e*. 

Let us explain more details. Our basic idea is to approximate PDE by discrete 
rational dynamics of the form: 

^N+l — J \^N-la ' • • • I -'AT I ^N-h ' • ■ • ' '^N+ki ) ■ • • i ^N+ka+i' 

by introducing scaling parameters 2;^ = t™'u{x, s) and {N, t) = (e^^x, e^'^.s). Notice 
that for any n variable function / as above, one can choose various types of the 
sets of variables {zll'\ , . . . , z%_i^ z%_^_^^ z*j^_^^^_^^ ) . We say that the 
sets: 
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arc the approximation data. 

Let us say that the above rational dynamics is consistent, if all fci , . . . , kd+i < 1 
hold. As a general procedure, the rational dynamics with the scaling parameters 
above give pairs of partial differential equations F(e, u,Ux, ■ ■ ■) = as the leading 
terms, and the error terms F^(e, u, Ux, . . . ) by use of Taylor expansions (3.S): 

+ e™~^^F^{e,U,Ux,Us,Uxs, ■ ■ • , Was, W(„+i)3;(^„+i,o)) • • • , W(a+i)s(^0,a+l)) 

We say that a partial differential equation P(m, Ux,Us, Uxs, • ■ • , Uax, u^s) is in PDE^, 
if there is an induced pair (F, F^) as above and a positive function /i > so that: 

(1) F satisfies the equality: 

r Ux, lis, Uxs, • ' ' •> Uaxi '^as/ — "7~7 T 

llyU, Ux, Us, Uxsi • • • 5 Ucjixi Uas ) 

(2) there is a constant C > so that the pointwise estimates hold: 

\F^{€,U,Ux,Ub, . . . , Uas,U(^a + l)x{^a+l,o),Uax,s{^a,l), ■ ■ ■ , U{a+l)si^O,a+l))\ 

a-\-l)x (?a+l,o)| + \Uax,s{^a,l)\ H H |M(a+l)s(^0,a+l) |)- 

P is said to be in PDE2*", if in addition the discrete dynamics is consistent. 
For the domains of solutions, we consider on: 

Uo,Ao) X [0,To) P e PDE^™, Ao,To € (0,oo], 

\(o,oo) X [o,ro) PePDE~, Toe (0,00]. 

The following result was obtained by the analysis of the combination of Lipschitz 
geometry with some scaling limits called tropical geometry: 

Theorem 1.1 (K2). Suppose two PDEs P,Q G PDEj are obtained from two 
relatively elementary and increasing functions f and g which are mutually tropically 
equivalent, where * = fin or 00. Then they are mutually related: 

P-^tQ i/P,gePDE^ 

where + = e or e^, and + can be e only when both f and g are linear. 
As a concrete case, we have the following: 

Proposition 1.2. Two partial differential equations of order 1; 

e 1 2 £ / 

Vs + -VVx--V =0, 2Us + -U(Us + Ux) = 

are both in PDE^ in the class (20, 2, 1, 2, 1), and they are mutually related in ~^ 
in the class (40,4, 1,2; 1). 

The following theorem suggests that "sizes of PDE2/ ~J will be large": 

Theorem 1.3. (1) For any M,D,L > 1, there are Iq and ao,bo so that for each 
even I = 2m > Iq, two PDEs u, v : (0, ao) x [0, bo) — >■ (0, 00) of order l: 

eus + —U2s H h j^uis + e^Uxs + —Uxss H h |yM(j-i)M's = 0, 

11 e' 
e{vs + Vx) + e^{Vxs + -^V2s + -^v^x) H h -j^Vix = 
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are both in PDEj " in the class (1,1,0,1,1), and they are mutually unrelated in 
''^fin in the classes {M,l,D,L;l). 

(2) For any M,c,D,L > 1, there is Iq so that for all I > Iq, pairs of PDEs 
given by: 

4us + eu{us + Ux) = 0, Vs= Ivx 

are both in PDE^ in (20, 2, 1, 7 + 1, 1), and are mutually unrelated in in the 

class {M,c,D,L: 1). 

The analysis for (2) touches with the technique in the field of the singular per- 
turbations in the sense of continuity of solutions at e = 0. 

The discrete dynamics used in the proof of (1) above satisfy the following prop- 
erty (def 4.1): 

Lemma 1.4. There is a pair of the discrete dynamics which are mutually infinitely 
unrelated. 

Now in particular one has obtained the foUowings: 
Corollary 1.5. There are unrelated pairs both in 

(1) (PDE{™, ~}„), (2) (PDE-, 
with respect to any constants. 

Finally we would like to address some problems which arise from the line of this 

paper: 

Question 1.1. Are there unrelated pairs in the following classes: 

(1) (PDE^», (2) (PDE,-, ~^), (3) (PDE,-, ~^). 



In order to approach such problems, computer calculations seem very convenient. 

In fact in many examples estimated here, one can express very concrete constants 
even though they are far from optimal. So in section 5, we have given computer cal- 
culation results. It turns out that certainly their data reflect related and unrelated 
situations. 

In order to treat the real models in physics, we need to extend this method 
for systems of partial differential equations. In particular the Korteweg-de Vries 
(KdV) equation is quite intriguing systems and has surprising relationship with the 
sohton cellular automaton system known as Box and ball system (BBS) through 
the time-discretized Lokta-Volterra system. We would expect that study of such 
direction might lead us to rough classiflcations of non-equilibrium systems. 



2. Classes for partial differential equations 

2. A Initial conditions: Let us consider C"+^ functions u, v : (0, Aq) x [0, Tq) 

(0,oo), where Ao,To € (0,oo]. 

For < e < 1, we introduce the initial rates: 

[U : V]e = sup {— r-)=^ . 

(x,s)e(0,Ao)x[0,e]U(0,e]x[0,To) ""{X, S) 
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Let u : (0, Aq) x [0, To) — )■ (0, oo) be a function of class C""*"^ as above. Then we 
introduce the uniform norm of u of order a + 1 by: 

Suppose u{x, s) > is positive. Then we introduce the higher derivative rates 
by: 

||u||a+l 



iiif(x,s)e(o,Ao)x[o,To) u{x,s) 



2.B Analytic relations on partial differential equations: Let: 

(*,+) e {(/m,e), (/m,e^), (oo,e), (oo,e^)}. 
Definition 2.1. Let P,Q £ PDE2 of order a. P and Q are related in {*, +): 

P Q 

if there are constants M,c,D,L > 1 and C,C' so that for any positive solutions 
u,v:iO,Ao)x [0,ro) ^ (0,oo).- 

P(e, U, Ux, Us,U:rs, • ■ • . Mas) = 0, Q{e, V, V^, Vs,Vxs, • • • , fas) = 

they satisfy the asymptotic estimates: 

v{x, s) - |m^"''(^+«+i) \u : v]Le + = e 

for all < Le < min(^,^o,^o,C''); where K = max(ifa+i(u), ifa+i(w)). 
The following is immediate: 

Lemma 2.1. Let P,Q € PDE^. Pi ~+ P2 holds whenever Pi P2 for both 
+ = e or e^. 

When we specify the constants, then we denote: 

Pr^fQ in {M,c,D,L; a), or P ~J Q in {M,l, D, L; a). 

The following is immediate: 

Lemma 2.2. Suppose M <M', c<c', D< D' and L < L' hold. Then Pi ~+ P2 
in (M, c, D, L; a) implies Pi ~+ P2 in (M', c', D', L'; a). 



Remark 2.3. Notice that even though two pairs Pi,P2 € PDEg are unrelated 
in ~+ in the class {M,c, D, L;a), they might still be related in ~+ in the class 
{M',c',D',L'-a). 

In [K2] , wc have obtained a general method to produce related partial differential 
equations. In the next section we sketch its construction, and in section 3.D.2, we 
calculate a concrete example of a pair of related PDEs which arise from the method. 
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3. Construction of related PDEs 

3. A Evolutional discrete dynamics: Let / = ^ be a rational function of n 
variables, where k and h are both polynomials. An evolutional discrete dynamics is 
given by flows of the form {z^}/,,Ar>0: where one regards that t is time parameter: 



^0 




• ■ 


. lit 


= 0) 


4 


4 ■ 


7I 


■ i{t 


= 1) 












4 




• -^N ■ 







A general equation of evolutional discrete dynamics is of the form: 

~ J \^N-lo' • ■ • ^ I ^Af-ii ) • • • ) ^AT+fci ) ■ • ■ ) -^N+kd+i ' 

where /j, kj > 0, N > max(/o, • • • , Id+i) and t > d, with initial values: 

^0 = {^a}o<a<max(io,.-->'<i+i)>*=0,lv ^ {^Jv}o</(<d,Ar=0,l,... • 

3.B Approximation of PDE by rational dynamics: Let us consider C"*"*"^ 
functions u : {0,Ao) x [OjTo) — )■ (0, oo) where Ao,To G (0,oo]. We consider its 
Taylor expansion, where \{x,s) — \ < \{ie^,je'^)\: 



-2jv+i = u{x + ie^, s + je«) =u + ie^u^ + je'^Us H ^"2x H ^ — 



+ je^iePw^, + • • • + 
(ieP)("+i) 



+ 



(a + 1)! 



(a + 1)! 

Let / = be as in 3. A and choose an approximation data. Then we consider 
the corresponding discrete dynamics z*j^^-^^ = f{z*j^^i^,. . . , J and change of 

variables e"^u{x, s) = z\^, N = and t = If < oo, then we assume fcj < 1 
for all 1 < i < d+ 1. 

Let us insert the expansion and take their difference: 

= e™w(x + eP,s + e«) - /(e'"M(x - loe^, s + e"), e"M(x + kd+ie^, s - de«)) 

_ c"'F\ti) + ('"-i'F-(ti,,) + ("'+''F-'((/,) -r (-"•+i'F'-{ti, a,,) + ... 
~ h{€'^u{x - loeP, s + e«), . . . , eu{x + kd+ieP, s - de")) 

where F^(u, u^, . . .) are monomials. 

Let us divide its numerator into two parts: 

e™Fi(u) + e"'+PF^{u^) + e"'+'^F^{u,) + €^"'+pf\u,u^) + ... 

a-\-l)x 

so that the first term contains only monomials whose derivatives of u are up to 
order a, and all the monomials of the second term contain derivatives of u of order 
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a + 1. Then putting: 

= P{e,u,Ux, . . . ,Uas)/h{e'^u{x - Iqc'', s + e^), . . .), 

= Q{e, u, .., {uixjs{^ij)}i+j=a+i)/h{€'^u{x - loe^, s + e"), ...), 
one obtains the following expression: 

J V^JV-/o ' ■ • ■ ' ^AT+fcd+i ) 
— 6 F(6, Its, Uxsi ' • * ? "^axt ^as) 

Definition 3.1. partial differential equation P{u, Ux,Us, u^s, ■ ■ ■ , Uax, Uas) is in 
PDE^, if there is an induced pair (F,F^) as above and a positive function h > 
so that: 

(1) F satisfies the equality: 

h[e"^u{x — loeP,s + ei),...) 

(2) there is a constant C >0 so that the pointwise estimates hold: 

\F^{e,U,Ux,Us, {Uixjs{^ij)}i+j=a + l)\ 

< C(|'U(a+l)x(Ca+l,o)| H 1- |W(a+l)s(^0,a+l)|)- 

P is in PDE2'"' if in addition, the discrete dynamics is consistent. 
Let us consider the last term with the estimate (see 2. A): 

C{\u^a+l)x{^a+l,o)\ H 1- |M(a+l)s(Co,Q+l)|) < CI ||u||a+l 

where I is the number of the summation of a + 1 derivatives of u. We say that 

the number CI is the error constant for the approximation of P. Notice that the 
error constant is determined by the approximation data and the original rational 
function /. 

For P e PDE2*", we consider (0,^o) x [0,To) as domains of solutions, and 
(0,00) X [0,To) for P e PDE^, where Ao,To e (0,oo]. 

Let: 

be the approximation data. Then one obtains several numbers: 

L = max(^,(i), D = max(p, g), A: = max(fci, . . . , /c^+i), ^ = max(Zo, • • • , ^d+i)- 

Below in 3.C, we induce extra numbers M = Mf and c = Cf from /. Then the 
function / and the approximation data determine the constants: 

{M,c,L,k,D). 

We say that P{u, Ux, Us,Uxs, ■ ■ ■ , Mas) is approximable in the class {M, c, L, k, D) in 
PDE2, if it is induced from some discrete dynamics as above whose constants are 
all less than or equal to M,c,L,k,D respectively. 
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Remark 3.1. Notice that three data, defining equations of discrete dynamics (1) 
z*j^^i = f{z*^^i^, .. ., -z^+fe^^j); (^) the exponents of the scaling change of variables 
{rn,,p,q) and (3) the order to take the Taylor expansions a, determine the defining 
PDEs. 

3.C Scale transform and tropical geometry: A relative (max, +)-function 

is a piecewisc linear function of the form: 

^p{x) = max(Q!i + dix, . . . , am + cimx) ~ max(/3i + bix, ...,/?; + bix) 

where aix = Yj^^-^a\xi, x = {xx, . . . ,Xn) G M", a; = (a^^, . . . , a"), 6 G Z" and 
ai,j3i e M. Notice that ip is Lipschitz, since it is piecewise linear. We say that the 
multiple integer M = ml is the number of the components of Lp. 

Correspondingly tropical geometry associates the parametrized rational function 
given by: 

' ~ ht{z) Y.i^^thzb'^ 

where z" = W^^-^zf , z = {z\, . . . ,Zn) € K>o- We say that ft above is a relatively 
elementary function. We say that both terms ht{z) = 'S^f.^-^t^'' z'"' and kt{z) = 
E^j^t"'=z'''= are just elementary functions. 

These two functions </? and ft admit one to one correspondence between their 
presentations. Moreover the defining equations are transformed by two steps, firstly 
by taking conjugates by logj and secondly by letting t ^ oo (sec [Mi] ) . 

In some cases the same (max, +) function admits different presentations, while 
the corresponding rational functions are mutually different. For example for ip{x) = 
max(0,2a;) = max(0, x, 2.x) = iIj{x), the corresponding rational functions ft{z) = 

+ 1 and gt{z) = + z + 1 are mutually different. 

Definition 3.2. Let ft and gt be the relatively elementary functions with respect 

to (f and ip respectively. 

We say that ft and gt are mutually tropically equivalent, if ip and ip are the same 
as maps (but possibly with the different presentations). 

For a relatively elementary function ft, let Cf > 1 he the Lipschitz constant and 
Mf be the number of the components with respect to the corresponding (max, +)- 
function (p. 

3.D Asymptotic estimates between different PDEs: Let / be a relatively 
elementary function of n variables. We have two constants Mf,Cf > 1 in 3.C. 
Let {{z'i+\,...,z'}+\z*^_i^,...,z%^k^,...,z*j^^^^^J,{m,p,q)} be an approxima- 
tion data, and L, k, D be the corresponding set of the numbers. With these num- 
bers, consider: 

M = max(M/, Mg), c = max{cf,Cg), C 

where C is the maximum of the error constants for the approximations of PDEs 
(3.-B). We have obtained the following result: 

Theorem 3.2 (K2). Let f and g be both relatively elementary and increasing 

functions ofn variables, which are mutually tropically equivalent. Let P, Q G PDEj 
be two PDEs of order a > which are approximated by f and g with the above 
approximation data respectively, where * = fin or oo. 
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Let us take positive C"+^ solutions u,v : (Q,Aq) x [0,To) — > (0,oo) with: 

P{e,U,Ua:,Us, ■ . .,Uax,Uas) = 0, Q{e,V,V3;,Vs, . . . ,Vax,Vas) = 

and assume the estimates Ka+i{u), Ko:+i{v) < K. Then for any < e < minl^^, {L+ 
l)~^Ao, (L + l)~^To,n~^), they satisfy the asymptotic estimates: 

// both f and g happen to be linear, then c = 1 and so they admit the exponential 
asymptotic estimates. 

The core of such asymptotic estimates has also appeared in [Kl] in the case of 
discrete dynamics. 

Corollary 3.3. Under the above situation, Pi ~+ P2 in (M', c', D, L+1; a), where: 

Um)^,c'^) 01, 
[((2M)8^l) c = i 

Proof: Firstly suppose c > 1. Then one has the estimates: 



(M',c') 



{2Mf — < [(2M)3= 



_ (c*^).— (x+,+1) 



Next suppose c = 1. Then: 

lim(2M)^^^ — = (2M)^(^"'°(^+'='*)+i) < [(2M)^*^]^"°(^+"+i). 

c— >1 

This completes the proof. 

3.D.2 Example: Let us apply the above method to a concrete case. See [K2] 
for more examples of order 2. 

Proposition 3.4. Two partial differential equations of order 1: 

Vs + ^VVx - = 0, 2Us + + Ux) = 

are both in PDE^ in the class (20, 2, 1, 2, 1), and they are mutually related in 
in the class (40, 4, 1, 2; 1) 

Proof: We have verified that these equations are in PDE^ in [K2]. For conve- 
nience, we give the explicit approximation data here. Let us choose: 

^t+l _ t+l t t ^1 = gjv(l + 2zy'_\) 

~ ./ ^JV) — 2 2(1 + ) ' 

= ff(w/_i, W;jv,W;v+2) = + 2(1 + ' 

and (m,p, q) = (1, 1, 1). Then it is direct to check the foUowings: 
ea{x + e, s + e) — A{ea{x — e,s + e), ea{x, s), ea{x + 2e, s)) 

_ f2(TTi^(2ws + 2ew^-t;2) + e2Fi(w,Wa:,",Uxx) {a,A) = {v,f), 
\2(^^{^'^s + euus + euux) + ^"^G^i^Ux, ..,Uxx) {a, A) = {u,g), 
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where the estimates hold: 

^ ^^"'^^1''^^^^-''^^ + e(2|z;,,(6i)| + k..(6o)l) < 5e||^;||2, 

, < en6Ks(r; n)| ^ ^ ^ |n..(r?2o)|) < 4e||n||2. 

(1 + eu) 

Thus both the equations are in PDE^, and their error constants are both bounded 
by 5e < 5. Thus combining with theorem 3.1, one has verified: 

Lemma 3.5. Let us put K = max.{K2{u) , K2{v)) and choose anyO <e< min(0.1/r~^, 1/3). 
Then any positive solutions to the above two PDEs satisfy the asymptotic esti- 
mates for all (x, s) € (0,oo) X [0,To); 

v{x, s) 

Proof of proposition 3.4; The above estimates are bounded by: 

for all < e < min(0.1ii'~^, 1/3). This completes the proof. 

4. Unrelated classes 

Let us take the discrete dynamics -z^i = /{z^nIi^,- ■ ■ , z^j^^^^^J and {m,p, q). 
Then one obtains a family of PDEs {Pa}a>i with respect to the order of the Taylor 
expansions (Remark 3.1 in 3.5). 

Let us fix {m,p, q), and take two discrete dynamics = f{z]^]^i^, . . . , -z^^fc^^^) 
a,n.dw*j^^i = g{z*N^ig, ■ ■ ■ > -2^J/+fc^_,.i)- Then one obtains two families of PDEs {Pa}a>i 
and {Qa}a>i correspondingly. 

Definition 4.1. Two discrete dynamics z'j^^^ = f{z^]^^i^^, ■ ■ ■ ,z^j^_^^^_^J andw*j^^-^ = 
g{z'j^^j^^, .... '-'/v+fe^^j) are infinitely unrelated in if for any constants M, c, D, L, 
there are some a so that Pa and Qa are mutually unrelated in ~+ in the class 
{M,c,D,L;a). 

Here we give some pairs of PDEs which are mutually unrelated. We treat two 
cases, where: 

(1) they are both in PDEj'", which are exponentially unrelated in They 

arose from two discrete dynamics which arc mutually infinitely unrelated. 

(2) they are both in PDE^, which are doublc-cxponcntially unrelated in 

4. A Unrelated pairs in the linear case: Below we show that the following 
linear equations are mutually unrelated: 

Theorem 4.1. For any M,D,L > 1, there are Iq and ao,6o so that for each even 
I = 2m> lo, two PDEs u, v : (0, oq) x [0, 60) (0, 00) of order I: 

eUs + —U2s H h -j^Uis + e Uxs + -^U^ss H h = 0, 

11 

e{Vs + Vx) + <^'^{Vxs + 2^2s + 2^2x) + 1- -jjVla; = 
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the fallowings; 

(1) they are approximable in the class (1, 1,0, 1, 1) in PDE2'", and 

(2) they are mutually unrelated in {M,l,D,L;l). 

Proof: (1) Let us consider the discrete dynamics given by: 



z 



N+1 ^ -^W+li "'at+I — '^N 



and the scahng parameters hy N = ^,t = j and z^f = u{x,s). Then by taking the 
Taylor expansions up to order l + l, one obtains the desired PDEs, which are both 
approximable in PDE2™ in the class (1, 1, 0, 1, 1, 1). Notice that each monomial of 
the PDE in u contains derivatives of s. 

(2) Let f{x) = + 1, and put u{x,s) = f{x) and v{x,s) = f{x — s). It is 
immediate to see that they arc solutions respectively, because of independence of 
the variable s for u, and of the symmetry of the equation for v. Moreover the 
equalities K{u)i+\ = Kij^i{v) = hold, since I + 1 derivatives of u and v are both 
equal to zero. We require that / are even in order to guarantee positivity of values 
of V. 

If they were related in (M, 1, Z), L; Z), then there is some C independent of solu- 
tions so that they must satisfy the asymptotic estimates: 

v{x, s) 

for all Q < Le < min(ao, ho,C). 

Let us choose small e > with Le< 1. Then the estimates hold: 

^e-°(x+s+l) ^ ^e-"(ao+6o + l)^ 

x' + l (Po-xY + l. 

[u : v\Le < [u ■■ v\i = maji{ sup , , , sup r— ) 

o<x<ao - Ij' + 1 o<x<i a;' + 1 

where: 

a;' + 1 
{x -iy + 1 '~ [1- 

{bo -xY + l , 

sup T— = 60 + 1. 

o<x<i a;' + 1 

Let us choose bo > 3. Then the estimate [u : v]Le < + ^ holds. 

Now let us find oq > &o ^ 3 so that the inequality ^^35^ > 60 holds. In fact for 
2 < /3 < 4 (say /? = 3 is enough), if one chooses large ao with Uq — /3ao > then 




bo = 2 («o + y - /3«o) > 3 

satisfies the required conditions. Notice that the equality = b holds, in the 

case when b = ^(a + vc? — 4a). Thus if one chooses bo as above, then the estimate 
> bo holds by an elementary observation. 



12 TSUYOSHI KATO AND SATOSHI TSUJIMOTO 

Let us fix such a pair (ao,foo)- Now = ^g_^'^^y_^-^ holds. If one chooses 

sufficiently large I » 1, then the estimates: 

u{aoM) .u 4 + 1 1 f^if (ao-6o<l) 

v{ao,bo)^' ' (ao-M' + 1(6^1) |> 1^^^ (ao - &o > 1) 

hold, since the last term is independent of I. Then one has: 
t;(ao,Oo) 

This is a contradiction. This completes the proof. 

4.B Non linear estimates: Let us treat the case of double-exponential esti- 
mates. Let P{e,u,Us,Ux, ■ ■ ■ ,Uas) = be in PDE^ of order a > 1, and let us 
compare its solutions with the translations: 

Pi:v, = Iv^ {I > 0). 

Let us start from a general situation. 

Lemma 4.2. Suppose that for some j > So > 0, C,C',C" > and for all small 

< e < eo; there are solutions u^{x,s) on {0,Sq) x [0,(5o) — > (0,oo) with the initial 
values u{x,0) = f{x) — 1~ x, which satisfy both the estimates: 

(1) C <u,{x,s)<C', (2) K^+i{u,)<C". 

Then for any M, c,D,L> 1, there is some Iq > so that for all I > Iq, P and Pi 
are unrelated in in the class (M, c, D, L; a) 

Proof: We verify the conclusion for a specific Iq > 0. Then the general case 
follows by restricting the domain of the solutions by choosing some smaller Sq > 
S',>0. 

Let us put v{x, s) = f{x + Iqs). Then v is the solution to Pjp and Ka+i{v) = 
holds. 

Let us take sufficiently small (5o>>e>>(5>0, and choose Iq with /o<5o = 

1 — S() — S. One may assume the estimate IqLc < Sq. Then Sq + loLe < 2So and 
Le + IqSo = l-6o-5 + Le <1- ^. So the estimate x + IqS < 1 - ^ holds on the 
initial domain {x,s) G (0,(5o) x [0,Le) U (0,Le) x [0,5o)- Thus the estimate holds: 

[v : u,]Le < 2C'6^\ 

Suppose they could be (M, c, D, L; a) related. Then they must satisfy the as- 
ymptotic estimates (^^^)^^ < M''^'''' {2C'Sof'' . 

On the other hand v{So, 6o) = S and u,{5o, Sq) > C hold. So > f which 

can be arbitrarily large. This is a contradiction. This completes the proof. 

4.B.2 Conservation equations: Let us apply the above situation to the first 

order conservation equations and the translations: 

AUs + €u{Us +Ux) , Vs = IVx- 

We verify that both are in PDE^ and are mutually unrelated in (Compare 
this with proposition 3.3): 
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Theorem 4.3. For any M, c, D, L, there is Iq so that for all I > lo, pairs of PDEs 
given by: 

4us + eu{us + Ux) = 0, Vs ^ Ivx 
are both in PDE^ in (20,2,1,/+ 1,1), and are unrelated in in the class 

{M,c,D,L;l). 

Before proceeding, let us briefly recall a way to produce solutions, called the 
method of characteristics for the conservative non linear equations of the form 

Us + F{u)ux = 0. 

Let 0) = /(^) be the initial condition, and try to solve the equation ^ = 
F{u) for u = u{x, s) with x{^, 0) = C Then ^^^^ = + = holds. So u is 
constant along x{^, s). Moreover = || = o, and so s) = F{f{^))s + C 
Thus if one could solve ^ = ^{s, x), then f{^{x, s)) will give us solutions, since u is 
constant along x{^,s), 

Proof: For the second equation, it is induced from the discrete dynamics w'pj'^-^ = 
so it lies in the class (1, 1,0,/+ 1, 1). Thus both are in PDE^ in the 
class (20, 2,1,/+ 1, 1) by proposition 3.3. 

Now let us consider the initial function / : (0, 1) — > (0, oo) by /(x) — 1 — x. The 
result follows if one can find solutions u which satisfy two conditions (1) and (2) 
in lemma 4.2. We will solve the equation very concretely by use of the method of 
characteristics. 

Let us consider the equation: 

^"i + M./(e)"i + Mi-e)' 

Then it gives the equation /x^-^ — (1 + l^ix — s + 1))^ + (1 + — /US = 0, and one 
can solve it with x{^, 0) = ^: 
1 



^{x, s)^— [ 1 + ^(a; - s + 1) - + M(a; - s + 1))^ - 4/x{(l + ^.)x - ^is} ]. 

Let us put the solution u{x,s) = 1 — ^(x, s). Then for a small (5o > and all 
sufficiently small e > 0, two conditions (1) C < u{x,s) < C and (2) K2{u) < C" 
are certainly satisfied, by elementary calculations. 

This completes the proof 



Remark 4.4. One has an observation of positivity of solutions, which rely heavily 
on the structure of the equation. Let us rewrite the equation as Us = — euwx(4 + 
eu)~^. At s = 0, Ux{x,0) = — 1 < hold and so Us{x,0) = eu(4 + €u)~^ > 
holds. Thus there is some Tq > so that Us{x,s) > still hold for all {x,s) G 
(0, i) X [0, To). In particular u{x, s) > ^ for all (x, s) € (0, i) x [0, Tq). The same 
argument works for the equation Vs + ^vvx — ^v^ = 0. 

5. Computational aspects 

Our mathematical framework is intimately familiar with the computer systems. 
Many essential quantities of the analysis are computable within finite steps and 
finite values. By use of numerical calculations, we reprove Theorem 1.3 for some 
particular values of constants (M = 10'^ and I = 10^). This can be done pos- 
sible since the discretization of the PDEs are the rigorous solutions rather than 
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just approximating. However in more general situations, such verification would 
be expected to work effectively by use of numerical simulation with guaranteed 
accuracy. 

For the proof that the two partial differential equations are of unrelated pair, it 
is sufficient to show that there exists the test point does not satisfy the inequality 
of the asymptotic estimates in the definition 2.1. We focus on this inequality and 
give some discussions from the point of view of the numerical calculations in two 
cases: the related pair and the unrelated pair. 

Recall that we have induced PDEs from discrete dynamics, and our computation 
here will be done for these discrete dynamics. 

Let {z%}N,t>o and {w'^}jv,t>o be two discrete dynamics. Then for L, No and to, 
we put the discrete version of the initial rates by: 

[{z%} : ^ sup max{(^)±, (^)±} 

Let us denote u{eN,et) = z]^ and v{eN,et) = w% respectively. Then we regard 
that both u and v approximate u and v respectively: 

u{eN, et) ~ u{eN, et) 

and similar for v. In fact in the examples in 5. A, both u and v coincide with u and 
V respectively, which induces Proposition 5.1 below. 

Throughout this section, we choose the rescaling parameters {m,p, q) = (1, 1, 1). 
All of the numerical calculations in this section are performed by using the com- 
puter algebra system "Maple 13" with rational or floating number manipulations. 
Then numbers after calculation are converted to the floating-point numbers for 
presentation purposes. 

5. A Case of the unrelated pair: Let us recall two PDEs: 

6^ 6*^ 
eUs + —U2s H h Ij-WJs + e^Wxs + yMiss H h |yU(i_i)xs = 0, (*i) 

11 

e(ws + "x) + (^^{Vxs + -^V2s + 2"2x) H H -j^Vi^ = (*2) 

which are mutually induced by the discrete dynamics: 

Recall that in 4. A we have verified that the pair of the solutions u{x, s) = + 1 
and v{x, s) = {x — s)' + 1 breaks the exponential bounds: 

, U(X,S) v(x,s) . , „-J3r-r-l-<.-l-l'lr 1 

v(x,s) u(x,s) 

Now for N,t = 0,1,2, ... , let us put the solutions to the above discrete dynamics: 

= {eNy + 1 , u;*^ = e'(iV - t)' + 1 

which are precisely the same as the solutions to the PDEs u{x, s) and v{x, s) by 
the scaling change of variables x = Ne and s = te respectively. 
We put: 



Qe{x,s) = logmax{^, ^} - log [{zj^} : {t«^}]l,,e-ix,e-ia 
whose values are much bigger than e~^{x + s + 1) logM in the above case. 
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We calculate the values Qe{x, s) at the points {x, s) — {Ne, te) for {N, t e Z | 1 < 
N < Ao/e, l<t< To/e}. 
Let us put: 

Qe{x, s) = max ^0, Qe{x, s)^ . 

Now we verify the following by computer calculations. These calculations are 
performed by using rational numbers which means that the calculated values are 
exact ones. 

Proposition 5.1. The pair oj the above equations (*i,*2) o-i'^ mutually unrelated 
m -^.„ in (10^, 1, 1, 1; lO^) at e = 1/2. 

Proof: Since the values of discrete dynamics and w]^ are precisely the same as 
the solutions z*^ = u{Ne,te) and = {x,s) = {Ne.te) respectively, it is enough 
to verify that Qe{x,s) certainly hit bigger values than e~^(x + s + l)logM = 
6(x + s + 1) log 10 at some points. 

Let us choose constants by / = 1000, L = and e = 1/2, and consider the values 
of Qe at each point [x^s) = {Ne,te) ior 1 < N,t < 8. For example (5e(l-5,2) = 

4 4 

logmax{^, ||} - log [{z%} : {w^}]o,3,4 

Then we show the values of Qe{x, s) in 4-digits of precision at the point {x, s) — 
{Ne, te) ioT {N,t < N,t <8^ 4/e} in Table 1. 

Table 1. Computed values of Q{Ne,te) (Unrelated pair) 



N\t 





1 


2 


3 


4 


5 


6 


7 


8 





0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


2 


0.0 


0.6931 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


3 


0.0 


404.8 


404.8 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


4 


0.0 


287.7 


691.8 


287.6 


0.0 


0.0 


0.0 


0.0 


0.0 


5 


0.0 


223.1 


510.1 


510.1 


223.2 


0.0 


0.0 


0.0 


0.0 


6 


0.0 


182.3 


404.8 


287.6 


404.9 


182.7 


0.0 


0.0 


0.0 


7 


0.0 


154.2 


335.8 


154.1 


154.2 


335.7 


154.0 


0.0 


0.0 


8 


0.0 


135.5 


287.0 


64.50 


0.0 


64.50 


287.0 


133.0 


0.0 



Now for D = 1,M = 1000, logM^ "i^+s+i) < \ogM' °(4+4+i) _ ^^24.3. So one 
can find that the exponential asymptotic estimates: 

Qe{x,s) < logM<^"°(^+''+i) 

does not hold at the several points. This completes the proof. 

Remark 5.2. Owing to the exact solutions of the partial differential equations, 
we are able to prove proposition 15.11 by using computer system. Nevertheless, in 
any case, even that we know little of solutions, it might be possible to apply the 
self-validating numerical method [Loh,0] in mathematical proof. 



5.B Case of the related pair: In case of the examples of the related pairs we 
present here, the numerical calculations do not directly provide the mathematical 
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proof of the relevancy. Nevertheless such calculations give us several insights on the 
actual behaviour for further analysis by comparison with the case of the unrelated 
pair. 

Let us recall two PDEs: 

Vs + ^VVx - = 0, 2Us + ^u{Us +Ux) =0 

which are mutually induced from the discrete dynamics: 



t+l _ ^Af(l + 

2JV+1- 2 + 2(1 + ^^) ' ^ ^' 



2{l+w%) 



+ 0/1 , ...t N (**2) 



In Proposition 1.2, we have seen the asymptotic estimates for all solutions, which 
are equivalent to: 

log ( _ 2^-'(-+2«)+3log([^, : v]2e) < 2^"' (^+2^)+^ log 40 
\v{x,s)J 

Now let us consider the solutions to the above discrete dynamics (**) with the 
initial and boundary values: 

z% = {eNY + 1 , w%^e\N-ty + l 

respectively, where {N, = {0} x N U N x {0}. 

Let u and v be as in 5. A, and u,v : [0, Aq] x [0, Tq] be solutions to the cor- 
responding PDEs respectively. Then as before we regard that both u and v ap- 
proximate u and V respectively at the points {x,s) = {NeP,te'^) = {Ne,te) for 
{N,tGZ\l<N <Ao/e,l<t< To/e }. 

We put: 

QAx,s) = logmax{4^, } - 2^-^(-+Wlog [{z%} : {<}]i,e-x,e-.- 
Here we calculate the values using floating numbers with 100 digits of precision. 



Qe' {x, s) = max ^0, (5e= (x, s) 



where we choose constants I = 1000, X = 1, e = 1/10 and Aq = Tq = 1. In 
particular the domain of PDEs arc {(a;, s) = (A^e, te) : < N,t < 10}. 

The numerical calculations verify that all the entries are equal to 0. In particular 
the estimates: 

Qeeix,s) < 2^+(^+2^)/nog40 
follows, which gives the numerical verification of proposition 3.3 for particular con- 
stants. Notice that approximately the value is given: 

2^+(^+'^)/nog40|,=i/io,.=.=ioe ~ 6.337 x lO'' 

So far we have checked that certainly double exponential estimates hold for 
these pairs. Next let us examine whether they might still satisfy the exponential 
estimates. 

Let us consider the values of Qdx, s) = max ^0, Qe^ {x, s)j , where: 
Qe{x, s) = logmax{^, ^} - log [{z%} : {w^5v}]i,e-ia;,e-is 
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with the same constants, I — 1000, L = 1, e = 1/10 and Aq = Tq — 1. Table 2 gives 
the result of numerical calculations. 

Table 2. Computed values of Qe{Ne,te) (Related pair) 



N\t 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 





0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


2 


0.0 


0.0 


0.095 


0.189 


0.258 


0.331 


0.393 


0.463 


0.524 


0.0 


0.0 


3 


0.0 


0.0 


0.196 


0.327 


0.461 


0.567 


0.682 


0.776 


0.0 


0.0 


0.0 


4 


0.0 


0.0 


0.229 


0.450 


0.618 


0.794 


0.926 


0.0 


0.0 


0.0 


0.0 


5 


0.0 


0.0 


0.318 


0.565 


0.813 


0.991 


0.0 


0.0 


0.0 


0.0 


0.0 


6 


0.0 


0.0 


0.334 


0.670 


0.899 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


7 


0.0 


0.0 


0.399 


0.677 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


8 


0.0 


0.0 


0.294 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


9 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


10 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



Let us compare their values with: 

101og40 = e-Mog40|e=i/io < e-\x + s + l)\ogAO 

< e-i(x + s + 1) log40|,=i/io,.=s=ioe = 30 log 40 

where the left and right hand sides are approximately 36.89 and 110.7 respectively. 
Thus the inequality above holds for any points {N,t) G [0, 10] x [0, 10]. 
Finally we would like to raise a question: 

Question 5.1: Are two PDEs: 

Vs + -Wx - 2^ =0' + -u{Us +Ux) = 

exponentially related in or in ^^j^ ? 
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